Diversity-Multiplexing Tradeoff in the 
Multiaccess Relay Channel with Finite Block 

Length 

Chung-Pi Lee and Hsuan-Jung Su 
Graduate Institute of Communication Engineering 
Department of Electrical Engineering 
National Taiwan University, Taipei, Taiwan 
Email: D96942016@ntu.edu.tw, hjsu@ cc.ee. ntu.edu.tw 

Abstract 

The Dynamic Decode-and-Forward (DDF) protocol and the Hybrid DDF and Amplified-and-Forward 
(HDAF) protocol for the multiple-access relay channel (MARC) with quasi static fading are evaluated 
using the Zheng-Tse diversity-multiplexing tradeoff (DMT). We assume that there are two users, one 
half-duplex relay, and a common destination, each equipped with single antenna. For the Rayleigh 
fading, the DDF protocol is well known and has been analyzed in terms of the DMT with infinite 
block length. By carefully dealing with properties specific to finite block length, we characterize the 
finite block length DMT which takes into account the fact that the event of decoding error at the relay 
causes the degradation in error performance when the block length is finite. Furthermore, we consider 
the situation where the destination does not have a priori knowledge of the relay decision time at which 
the relay switches from listening to transmitting. By introducing a decision rejection criterion such that 
the relay forwards message only when its decision is reliable, and the generalized likelihood ratio test 
(GLRT) rule at the destination that jointly decodes the relay decision time and the information message, 
our analysis show that the optimal DMT is achievable as if there is no decoding error at the relay 
and the relay decision time is known at the destination. Therefore, infinite block length and additional 
overhead for communicating the decision time are not needed for the DDF to achieve the optimal DMT. 
To further improve the DMT, we propose the HDAF protocol which take advantages of both the DDF 
and the Amplified-and-Forward protocols by judiciously choosing which protocol to use. Our result 
shows that the HDAF protocol outperforms the original DDF in the DMT perspective. Finally, a variant 
of the HDAF protocol with lower implementation complexity without sacrificing the DMT performance 
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is devised. 

I. Introduction 

In recent years, cooperative communication has received significant interest as a means of 
providing spatial diversity when time, frequency, antenna diversity are unavailable due to delay, 
bandwidth and terminal size constraints, respectively. Cooperative techniques provide diversity 
by enabling users to utilize one another's resources and has been extensively studied for the single 
source from outage probability analysis or diversity-multiplexing tradeoff (DMT) perspective [HI, 

Practical communication systems usually involve more than one user. One of the most typical 
models is the multiple-access channel (MAC). The capacity region of MAC is well known and 
the DMT is also developed in ||4] for MAC. In [[3||, B, cooperative diversity was extended to the 
multiple users cases. For most cooperative protocols, substantial coordination among the users 
are required, which may be impractical due to cost and complexity consideration. Alternatively, 
we consider the multiple-access relay channel (MARC),i.e. the MAC with a single shared relay 
and focus on the dynamic decode-and forward (DDF) protocol [3]. For DDF protocol, the relay 
does not decode until it is possible to successfully decode source information message. The 
relay then re-encodes the message and transmit it in the remaining coding interval. In this model 
we concern, the users need not be aware of the existence of the relay. All cost and complexity 
are placed in the relay and destination. Such an architecture may be suitable for infrastructure 
networks, where the relay and destination correspond the station having more resource(i.e. base 
station). Moreover, sice a single relay is shared by multiple users in the MARC, the extra cost of 
adding the relay per user may thus be more acceptable. Finally, for further enhancement of DMT, 
we propose a hybrid protocol which combines DDF and multiple-access amplify-and-forward 
(MAF) [|6||, [|7]| protocol to improve the diversity gain at high multiplexing gain region. 

A. Related Research 

The MARC was first introduced in H. In MARC, the relay helps multiple sources simul- 
taneously to reach a common destination. Information-theoretic treatment of the MARC has 
focused on two aspects, namely, the capacity region and the DMT (the outage behavior of slow 
fading channel in the high signal-to-noise (SNR) regime [9]). The achievable rate for the MARC 
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has been proposed in [fTOl . [fTTI . [fT2l . However, the capacity region of general MARC remains 
unknown. The DMT for the half-duplex MARC with single antenna nodes is studied in [6J, 
[Q, ESI, Il2l. In del, Q, it is shown that MAF protocol is DMT optimal for high multiplexing 
gains; however, this protocol remains to be suboptimal for low multiplexing gains compared with 
DDF strategy [|T3l . In addition, the extra overhead to communicate the channel realizations of 
source-relay links to destination is required for MAF. Another relaying strategy for the MARC is 
compress-forward (CF) [Q. In CF, the relay exploits Wyner-Ziv coding to compress its received 
signal and forward it to destination. It has been shown in [j2l that CF also acheves the optimal 
DMT for high multiplexing gains but suffers from diversity loss for low multiplexing gains, 
moreover, CF costs much larger complexity. 

The present paper focuses on the DDF protocol for the half-duplex, single relay, single antenna 
case due to its nice balance between complexity and good performance at low multiplexing gain. 
Moreover, we propose a hybrid DDF-MAF (HDAF) protocol to enhance the DMT by improving 
the poor performance of DDF at high multiplexing gain. 

B. Summary of Results 

Previous work, [fT3l . [|3l assume an infinite block length such that there is no decoding error 
at the relay and the number of relay decision moments is also assumed infinite. Inspired by [fT4l 
for the single user, single relay case, we analyze explicitly the achievable DMT with finite block 
length and finite relay decision moments for the MARC. Note the former is a special case of the 
MARC when there is only one user. Moreover, different from the proof of DMT achievability 
in lO, we do not separately average the probability of event at relay and destination over the 
random codebook since they are not fully independent, which will be made more clear in Section 
IB-C[ As [14J, two issues are discussed in our model: 1) the effect of decoding errors at the relay. 
2) the fact that the relay decision time is not generally known priori at designation. In order 
to tackle 1) a decision rejection criterion at the relay, such that the relay triggers transmission 
only when its decoding is reliable. For 2), the destination jointly decodes the relay decision 
time and the information messages based on the generalized likelihood ratio test (GLRT) rule. 
Our results show that in order to achieve the DMT, additional protocol overhead informing the 
destination about the relay decision time is not necessary and the loss of DMT due to decoding 
error at the relay can be avoided. Finally, we propose HDAF protocol which takes advantage of 



Fig. 1. MARC model 

both DDF and MAF protocol. Our analysis shows that with the finite block length MT, HDAF 
protocol outperforms the original DDF protocol, especially at high multiplexing gain region 
and low multiplexing region when M is small. HDAF also has better DMT than MAF protocol 
when M is moderate large. In addition, without causing any loss in DMT perspective, a variant of 
HDAF with lower complexity by allowing the relay's transmission only after M/2 relay decision 
moments [11511 is devised. Notice that the analysis with finite T will be much complicated since 
we have to deal with events which depends not only on channel realization(i.e outage) but also 
on codebooks when applying random coding techniques. 

In Section ini we introduce the system model and review relevant previous results. In Section 
Uni a characterization of the DMT of the DDF protocol for MARC with finite block length is 
presented. In Section |lVl we devise HDAF and its variant and characterize the DMT with finite 
block length. Section |V] gives the conclusion. 

II. Notations and Definitions 

A. System Model 

We Consider the two-user(Sl and S2) MARC model, a relay node (R) is assigned to assist 
the two multiple access users, (see Figil). The users are not allowed to help each other (due to 
practical limitations, for example). The relay node is constrained by the half-duplex assumption, 
i.e. the relay can not transmit and receive simultaneously. Each node equips with single antenna. 
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All wireless links are assumed to be frequency nonselective and block fading, where the 
channel coefficients are random but remain constant over the whole duration of a codeword and 
the channel coherence time is much larger than the allowed decoding delay. Let hi, gi, g,- denote 
the fading coefficients between the user i to relay, user i to destination, and relay to destination, 
respectively. The channel fading coefficients are i.i.d. C9{jS),l) variables, corresponding to i.i.d. 
Rayleigh fading. Here we assume the perfect channel state information at receiver (CSIR). We 
adopt the slotted transmission where a codeword spans M slots and each slot consists of length 
T symbols. There are thus a total block length of MT. In decode-and-forward protocols, the 
block of length MT symbols is split into two phases. In the first phase, the relay receives the 
signal from the source until the end of a certain slot, referred as the decision time. Then the 
relay tries to decode the source message. In the second phase, based on the decoded message, 
the relay sends its codeword to the destination in the remaining block . For DDF protocol, the 
decision time depends on the channel coefficient and the received signal. For the first phase, the 
signal received by the relay is 

2 

yr,k = Y,hiXi^k + nk, A:= l,2,...,mr (1) 

Note the decision time at the end of slot m is a random variable. The signal received by the 
destination is 

2 

yd,k = Y,SiXi,k + Vfe, ^ = 1 , 2, . . . , mr (2) 
During the second phase, the signal received by the destination is 

2 

yd,k = Y.8iXi,k + 8rXr,k + Wk, k = mT+l,mT + 2,...,MT (3) 

i 

We let Xi = [xj^i, ...X(^Mr]^ denote the user / codeword with rate Ri bits per symbol, is similarly 
defined but only the last (M — m)r symbols are transmitted. The same average power constraint 
P per symbol are imposed on each user and relay, 

where E[ ] denotes expectation over all codewords. The noise at the relay and destination are 
independent Gaussian noise with variance of a^, o^, denoted as rij^ ^ C1A^(0, a^) and ~ 
C!^C{0,ol) respectively, p = and p' = define the SNRs of the source-destination and the 
source-relay links ,respectively. In our model, we consider the symmetric case where Ri=R2 = 
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R/2, the sum rate is R bits per symbol. As [fT3l . the relay decision time is chosen to be the 
instant that the transmitted rate is within the achievable rate region at the end of m slot, which 
satisfies 

MTRi<m\og{\ + \hi\^p')T (4) 

MTR2<m\og{\ + \h2\^p')T (5) 

MTR < mlog (1 + |2 + \h2\^)p') T (6) 

m is set to the minimum m = 1,2,...,M— 1 such that ([5]), ^ hold, otherwise, m = M and 
the relay remains silent. 

For later use, we denote the complement of an event A by A, the transpose and Hermitian 
transpose of z by z^, z^, respectively. = x ii x> 0, otherwise equal to zero. Let x^,^, and 
x^„, denote the source transmit signal from time nT to kT and the relay transmit signal from 
time nT + 1 to kT respectively. y^„,, x^^, are similarly defined. 

B. Diversity-Multiplexing Tradeoff (DMT) 

Our work use a lot the notion of DMT posed in [|9]. We only provide definitions here. Consider 
a family of codes, such that the code has a rate of i?(p), corresponding to SNR p bits per channel 
use(BPCU) and error probability Pe{9) The multiplexing gain r and the diversity gain defined 
as 

,4lta5M, ,4-lto!flMP) (7) 

p^co log p p^oo log p 

and we can write as Pe{9) = P '^^ where = denotes the exponential equality . < and > are 
similarly defined. For the point to point multiple-input multiple-output(MIMO) channel with m 
transmit and n receive antennas, r < min(m,n), the optimal d*{r) is referred to as DMT. 

in. DMT OF THE DDF PROTOCOL FOR MARC VV^ITH FINITE BLOCK LENGTH 
A. Effect of Finite Block Length 

The DMT of DDF protocol for MARC with symmetric rate R\ =i?2 = f logp has been shown 



in mi. 



2-r j>r>0 



^DDF-MARC 



3(1 -r) i>r>i (8) 



3 ^ ' — 2 
2^ l>r>i 

r — — 3 
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However, to achieve the DMT in ([8]), a scheme with M oo possible decision times is necessary 
[[T3|. furthermore, an infinite block length T — t- oo is assumed such that there is no decoding error 
at the relay. For the practical code design, the code length is finite and the error event of decoding 
at the relay occurs even though the transmitted rate falls in the achievable rate region, i.e. (H), (|5]), 
^ are satisfied. Forwarding the wrong source information message would significantly degrade 
the error performance at the destination. Thus the probability of relay decoding error dominates 
the error probability at the destination. The DMT analysis with finite M and T has been treated 
in [fT4l for the relay case, a special case of MARC model if there is only one user. In our 
MARC model, there are two users interfering with each other, thus the derivation of DMT is 
more involved, moreover, we take an approach different from [[T4l to averaging error probability 
over all the random codebook. 

B. Characterization of DMT with Finite Block Length 

In this section, the proof uses the machinery of [O, [[T3l . H, jSl. Therefore, we only provide 
a sketch of the previous results involved and focus on the novel parts. First, we find an upper 
bound on the DMT by letting T oo, and assume the destination has the knowledge of m. 
By relaxing the constraint that T is finite, we characterize the DMT using outage probability 
analysis, hence make it as an upper bound for the finite block length. This is established by the 
following theorem. 

Theorem 1: The DMT of the two users, single-relay, single destination DDF scheme with 
decision times m = 1,2, ...,M and finite slot length T is upper-bounded by 

dout{r)= min {dmMf) + ^i^iAr)} (9) 

l<m<M 

where dm,R{r) is defined in (l40l) . and d,„,Dir) is defined in (|44l) - (|47l) . 

Proof: We use similar techniques developed in flU, [[T4l. Let Pout{r) denote the outage 
probability at the destination. P"l,t{r) denote the outage probability at the destination for a given 
m = m. The outage event for a given m = m is defined as 

(yE= {{gl.S2.gr) : (orz>U^2:DUO0,2),D)} 

where Oj"^, O^^, 2) d ^re defined as 

a^^D = {{gi.gr) ■mT\og{\ + \g,\'^) 

(11) 

+ (M - m)nog (1 + (1^1 12 + |g,|2) p) < MTR, ] 
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(12) 



^2D similarly defined by replacing 1 with 2. 

O'llalD ={igUg2,gr) ■ mriog (1 + + |^2|') p) 

+ (M - m)riog (1 + (1^1 12 + \g2\^ + p) < MTR} 
Then write 

M 

Pout{r)='£P{m = m)PZ{r) (13) 

m=l 

Since scaling SNR by a constant (p'/p) does not change the DMT. Define 

P^"^^^(r) = p-^'".i'W, l<m<M (14) 

/'(m = m) = p"^'»'«W, \<m<M 
Then Q clearly follows from (fT3T ) The proof of (fT4] ) is provided in Appendix |Al It remains to 
show that it is an upper bound of finite block length case. From standard arguments based on 
Fano's inequality it can be seen that Pgutif) is indeed the best we can get, thus completes 
the proof. ■ 

Note in the proof above, the relay decision time m is assumed priori known at the destination. 
The following theorem shows that this assumption is not necessary and the upper bound is 
achievable with finite T. 

Theorem 2: The upper bound of Theorem [T] is achievable for finite-length T and no priori 
knowledge of decision time at the destination decoder. 

Proof: In order to prove the achievability, we use standard random coding argument with 
bounded distance decoder [[T3l at the relay to overcome the effect of relay decoding error [14]. 
Owing to the introduction of the bounded distance decoder at the relay, the probability of relay's 
decision time at m slot and the probability of decoding error at destination given relay decision 
time m are not independent(the relay's decision time not only depends on the source-relay links 
but also on the codewords). Therefore, different from [[T4l where these two terms are separately 
averaged over the random ensemble, we take an approach to directly averaging the resulting 
error probability at the destination over the random ensemble. 

Codebook Generation: For given M, T, Ri = R2 = §logp, according to i.i.d components 
(riA^(0,/'), independently generate p^ codewords, Xi C C^^, for each / = 1,2 and X,- C 
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of cardinality p''^^. We let x,(w,), x^(w) denote the codewords in Xj and in Xr respectively, 
corresponding to the information message w,- G {l,...,p^}, w G {1, p''^'^}. 
Relay Decoding: From dH)-®, we define the relay outage event at slot m as 

o^={ihM) : (o;:^uo2:^u 0(1,2),^)} (15) 



and 



is similarly defined. 



OT,^^{h:\h\^<^^-J-} (16) 



Oil,2),R = \ (^i'^2) : \hi\^ + \h2\' < !> (17) 

Due to the finite block length MT, the relay may decode in error even (/ji,/j2) ^ O'^. Then 
an incorrect codeword is sent to the destination, causing significant interference. Thereofre, 
similar to [[T6l . [fT4l for the relay case, we introduce a bounded distance relay decoding decision 
function \\r^ defined as follows: for m= 1,...,M— 1, define the regions 5ra(wi,W2) consisting of 
all points y G C"'^ for which (wi,W2) is the unique message enclosed in a sphere of squared 
radius mr(l+5)a,^ centered at y, i.e., \y — hix'^Q{wi) — h2X2Q{w2)\^ < niT{l + 8)o^. If 

1) {hiM)iO^, 

2) y™Q G 5m(H'i,M>2), 

Then, Vs(y"})7^i7^2) outputs the decoded message w, where w = (vi>i,m>2) and the relay starts to 
transmit the signal x^^(m>) for the remaining part of the block. Otherwise, it waits for the next 
slot. 

Destination Decoding : The destination is not aware of the relay decision time m, hence it 
simultaneously detects the decision time and the information message according to the GLRT 
rule: 

{w,m} = argmaxp(y^o|w,m) (18) 

w,m ' 

where p{y^Q\w,m) is the decoder likelihood function. 

Error Probability Analysis : Let E denote the decoding error event at the destination and Er 
denote the decoding error event at the relay. Follow the steps in [[I4|. we have the following 
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Fig. 2. The DMT of the DDF protocol for MARC with finite M 



results, 

M 

= /'(m = m)P(£|in = m) (19) 

m=l 
M 

< P{jn. = m) {P{Er\m = m) + P{E\Er,m. = m)) (20) 

For m = 1, ...,M— 1, let 5 = ^ulogp, we have 

P{Er\m = m)<p-"''^^ (21) 

Since P(E\Er,m = m)>^^'^"''^^^\ we can choose a sufficiently large finite /j to make the terms 
P{Er\va = m) exponentially irrelevant in (|20] ). i.e. since dm,D{r) < 3, we can choose iIT > 3. As 
for the other terms in (|20] ). as mentioned earlier, they are not independent event since the relay's 
decision time not only depends on the source-relay links but also on the codewords. In Appendix 
IHl unlike [fT4ll . by averaging these two terms together, we show the average probability of the 
term P{m = m)P{Er\'fn = m) using Gaussian random ensemble is exponentially upper bounded 
by dout{r), thus completes the proof. 

■ 

It appears intractable to obtain a closed form of dD. In Figj2l doui{r) are plotted for M = 
2,5, , 10,20. As M grows, dout{r) is seen to approach the optimal DMT in [13] where M — )• oo. 
FiglHalso shows that even for a moderate value of M, doiit{r) achieves the DMT close to optimal 
one. 
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IV. Hybrid Amplified-forward and decode-forward protocol 
A. The DMT of Hybrid Amplified-Forward and Decode-Forward protocol 
It has been reported in [|6l [|71, the DMT of the MAF protocol is given by 

[2-f, 0<r<f 
dMAF{r) = < (22) 
[3(l-r), f<r<l 

The optimal diversity gain for high multiplexing gain (|<r<l)is achieved by MAF protocol. 
On the contrast, for the DDF strategy, the relay will only be able to during a small fraction 
of time slots, hence suffers form a loss compared to MAF. Combined with MAF protocol, the 
DMT may be further improved. This motivates us to propose a new hybrid strategy as follows: 

(1) For I < r < 1 , the relay simply uses the MAF protocol. 

(2) For 2 < r < |, the relay simply uses the DDF protocol. 

(3) For < r < J, the relay dynamically decode the source messages before M/2 (including 
M/2) time slots ,where we assume M is even. If it can not successfully decode , then the MAF 
protocol is instead used after the time slot M/2. 

The goal of the hybrid strategy aims to take advantages of both MAF and DDF protocol, 
we refer it as HDAF protocol. Note HDAF has the same setting as Section UlI] (i.e. finite block 
length) when DDF is chosen to be used. We have the following theorem and Figj2] showing that 
HDAF outperforms DDF especially at | < r < 1 and < r < j when the number of time slots, 
namely M is small since HDAF can achieve the optimal DMT of DDF protocol with infinite 
M for < r < J even a finite M time slots being used. Furthermore, HDAF has better diversity 
gain than MAF for ^ < r < | when M is large enough since it will approach close to the DMT 
of DDF with infinite M. 

Theorem 3: The DMT of HDAF protocol is 

2-r 0<r<^ 



doutir) ^<r<| (23) 



2^-3 
3(l-r) f<r<l 

Proof: To characterize the DMT of HDAF protocol dnoAF (r), for f < r < 1, the MAF 
protocol will be used and dnDAFi^) = dMAF{r). The advantages of MAF for high multiplexing 
gain is preserved. For 5 < r < |, the DDF protocol will be used and dHDApir) = dout{r) in 
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hence this region entails no loss in terms of DMT compared to DDF protocol. For <r<j, from 



[1131 . we know 2 — r is indeed the upper bound for MARC, we will show that dnDAFif) > 2 — r, 
thus establish duDAFi^) =2 — r. First consider the case where source message is decoded at the 
relay before M/2 time slots. In this case, since the upper bound of ^ is achievable by finite 
code length, we can similarly write dnDAFir) as 

dHDAFir) = ^ min {d!^,T^{r) +<^^^(r)| (24) 

l<m<M/2 I- 'J 

It clear that dj^^^^{r) = d,n,R{r)^ '^m.D^^(^) = d,n,D{r) for 1 < m < M/2. From (gD, the diversity 
gain at least 2 — r is obtained since dj^^^{r) = 2 — r. However, this is not necessarily the case 
since the source message may be decoded only after M/2 time slots and the MAF protocol will 
be used instead. It remains to derive the DMT of HDAF protocol for the case where decoding 
moment m > M/2. In [[6]|, the error provability Pe using MAF protocol is shown to be upper- 
bounded by 

Pe<Pei+Pe2+%.2) (25) 

where Pgj, 1= 1,2,(1,2) represent the probability of the error event that user(s) / are detected 
in error. Pg/s are averaged over the ensemble of channel realizations, and thus leads to lower 
diversity gain 2 — y for < r < |. However, for our case, the MAF protocol is used only when 
the relay can not decode the source message before M/2 time slots and the resulting diversity 
gain is expected to be larger. To deal with the HDAF case, from (l48l) . the probability that the 
source message can not be decoded before M/2 time slots can be split into two events 

m')( mm 

Then we need to upper bound the probability that MAF is used and decoding error occurs at 
the destination, defined as 

M 

|2 



PEMAF=P\EA{hxM) e Oj 



(26) 

M M ^ ^ 



^pyE,>^{{hxM)ioi}Ayl^iUM\ 

For the infinite block code length case, only the first term on the right hand side (|26] ) need to 
be considered. Nevertheless, in our finite block code length case , the second term need to be 
taken into account which much complicates the analysis. We upper bound the first term here 
and leave the analysis of the second term in Appendix O 
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Similar to (|25] ). partition the first term on the right hand side of (|26l ) into three mutually 
exclusive error events. To ease the notation, we denote them as Pe, in this Section. Define 
\hi\^ = p^"', \gi\^ = p^^', where i = l,2,r. From [3], the Pg^ conditioned on hi, gi is 

and define the outage event as 

(max{2(l-|3i),l-(|3, + ai)})+<r} 



(28) 



From ©-([n]), 

M 

o 



J=|a,:(l-aO+<-| /=1,2 
0(^12),^?= {(«bOC2) : (1 -min(ai,a2))+ < r} 
Let Pe, = P ''"i^'^', by discussing cases of (|28l) for different events of , it can be shown a 



diversity gain at least 2 — r is obtained. 



d.{r)= inf {ai+a2 + |3i+|32 + pj>2-r (29) 

dp^ (r) can be derived in identical manner. For dp^ (r), the DMT exponent conditioned on 
/ii,/z2 [,6j, is 

2(1 -r)+ min{ai,a2} > (l-r)+ (30) 

[3(1 - r) -min{ai,a2}]^ < min{ai,a2} < (1 
Then 

\y.2)= ^{«i+«2 + ^p,^^^j|(/zi,/z2)(r)}>(2-r) (31) 

(ai,a2)e07 

■ 

Theorem |3] shows that even for small M, HDAF can still achieve the optimal DMT of DDF 
protocol with infinite M when < r < ^. The loss in DMT of finite M compared to infinite M 
case comes from the fact that the event of erroneous decoding at the destination when the relay 
starts transmission after M/2 time slots becomes the dominant error event, which will not occur 
in the HDAF strategy since MAF is used instead in that case. 
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B. Variant of HDAF protocol 

Indicated in [fTSl . allowing the relay node to start transmission at any time slots may result 
in higher complexity since this requires a very high-dimensional constellation to ensure the 
possibility of source message being uniquely decodable within a small code length. We can also 
prove the modified HDAF protocol that allowing the relay to transmit only after the M/2 time 
slots do not entail any loss in the DMT perspective. 

Theorem 4: The modified HDAF protocol still achieves the DMT of Theorem [3l 

Proof: Note DDF protocol is used only for < r < |. Replace (l48l) with 1 to upper bound 
P{m = m) for all 1 < m < M/2, i.e. relay decodes before M/2 time slots. Following the remaining 
steps thereafter, we can obtain an upper bound dm,D{r). Observe (|44]) -(l46l). dm^oir) is already 
equal to optimal value 2 — r, 3(1— r)for0</<2,0</< 2^0rj' respectively. Since 2{3r^-^i) — 
5 for I > r > ^, decoding at / > ^ is good enough, we conclude that the modification does not 
affect the DMT achieved by the HDAF protocol. ■ 

V. Conclusion 

In this paper, we consider the design of cooperative protocol for MARC consisting half-duplex 
nodes. We analyze the DMT of DDF and proposed HDAD for MARC model with finite block 
length. The analysis captures the practical issue that the relay may decode erroneously even when 
the channel is not outage and the destination have no priori knowledge of the relay's decision 
moment. Our results show that additional protocol overhead informing the destination about the 
relay decision time is not necessary and the loss in DMT due to decoding error at the relay 
can be overcome. The difficulty of analysis comes from that we need to properly manipulate 
the events related to not only channel realization(i.e outage) but also codewords when applying 
random coding schemes with finite block length. We also propose HDAF protocol whcich shares 
both the advantages of DDF at low-medium multiplexing gain and MAF at high multiplexing 
gain region. It achieves higher(optimal) diversity gain than DDF protocol, particularly at high 
multiplexing gain and at low multiplexing gain when the number of time slots M is small. 
HDAF also outperforms MAF at low-medium multiplexing gain region when M is reasonable 
large enough. Finally, we have investigated the variant of HDAF with reduced complexity by 
allowing the relay to switch to transmission mode only after half of codeword and prove there 
is no entailing loss in terms of DMT compared to the original HDAF. 
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Appendix A 
Proof of Theorem [H 

Since d,n,R{r) is solely a function of the source-relay links and d,n,D{r) is a function of relay- 
destination and source-destination links. These two terms can be analyzed separately. 

A. Analysis of dm,R{r) 

We partition the event {m = m} into the set of events A"\ i.e.{m = m} = [jA"\ where / 
denotes any nonempty subset of {1,2}. For notation convenience, we make two extra definitions 



that are only used in Appendix lA] and Appendix |Bl Ofj^ = (|), = (|) where (|) denotes the empty 
set. Then for m = 1,...,M, A'p are defined as follows 

AT = {{huh2):0-^'PlO^} (32) 

Event A™ is similarly defined with 1 replaced with 2. 

^ri,2)={(^i'^2):o-72;,^n^} (33) 

Af denotes the event that source message can be decoded at slot m and the outage event O™^^ 
occurs. Note 

PiAT*) < P{m = m)= PiljAT) < ^/'(Af ) = P{Af.) (34) 



where /* = argmin lim — \ ' . Thus P{m = m) = P{A'Jl). Since | are exponentially dis- 

/ p'— >oo °" 

tributed and p = p', we compute P{Af) = p'^-*/"^'^'' as follows: 

1 ) Computation of dA'n (r), for M —I >m>2 : For r >0, M — 1 >m>2 



P{A'l') = { (35) 

Mr 



Due to the symmetry, dA>:^{r) = dA"^{r) is clear. 
2) Computation of dA"^^{r), for M — \>m>2: 



f >r>^, dAi^^ir)=0 (36) 
^>r, J^»,^(r) = 2(l-f^) 



m— 1 
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3) Computation of dm=i,R{r): 



P{m=\) 



m 
M 



> r, dm^x^R{r) =0 



4) Computation of dm=M,Rir): 



P(Af) 



^^^^ >r, J^M(r) = 1- 



dj^M{r) = d^M{r) is clear. 



M \ 
1,2) 



{M-i) 



M 



>r, J^M(r) = 2(l-Ji^) 



Collecting the results (I35l) -(|39]). overall, we obtain 



(37) 



(38) 



(39) 



dm,R{r) 



2(m-l) J _ 1 Mr 



0< r< 



3M 



3M 

<r< 



M ' 



2(m-l) 

J,„^^(r)=2(l 



OT— 1) 



^ ^ ' ^ ^7 a,n,R\' ) — ^ 



^<r<l, J„,i^(r)=oo 



(40) 



B. Analysis of dm,D{r) 

Similar to (|34] ). P"ut{r) = P(yn ^{r) the DMT d(yn^{r) has been derived in [fT3l . we provide 
results here, let / = 



and 

for i > r > 0, 



i>/>0, 2-r 



1 - - 

^ 2 



>/>i 2 



2(1-/) 



(41) 



!>/>!- 



2' 2/ 



<in»> (r) 

{\a.),D^ ' 



§>/>0, 3(l-r) 
l-r>/>f, 3-1^ 



(42) 
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for 1 > r > i. 



^>/>0, 3(l-r) 

^^.)z,W=r ~ (43) 
'''''' '!>/>§, 2i-i: 



Note dm,D{r) = min< d(yn {r),d(yn {r),dan (r) \. By some manipulation, we obtain. 



for i > r > 0, 



for if > r > i. 



forf >r> If, 



dmfi{r) = < 



i>/>0, 2-r 

i 2 ^ / - 2' ^ 2(1-/) 
2' TT 



1 > f > 1 _ r 2^ 



dm,D{r) 



2g^>/>0, 3(l-r) 



>/- 2(3r-l) ' ^~ 



l>/>l-i 



2-r 
2/ 



(44) 



(45) 



5r-2 



2(3r-l) 



>/>0, 3(l-r) 



2- 



5r 



dm,D{r) = < 



2-1+ 



>/> 
>/> 



5r-2 



2(3r-l)' ^ 2(1-/) 

2(l-r) 
2 ' / 



^ 4 



l-5>/> 



2- 



2(1-/) 



!>/>!-§ 



2-r 
2' 2/ 



for 1 > r > |, 



dm,D{r) 



f>/>0, 3(l-r) 
!>/>§, ^ 



(46) 



(47) 
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Appendix B 
PROOF OF THEOREM [2] 

We now turn to the analysis of the term P{m = m)P{Er\im = m) in (|20|) . First, we upper bound 
these two terms by extending the techniques in [|14I to our MARC model. Then , we average 
the product of these two terms over the Gaussian random ensemble 

A. Analysis of P{m = m) 

Let Zlin = Uvv=i Sm{w), where w = (wi,W2)- Averaged with respect to the random coding 
ensemble, we may choose without loss of generality wi = 1,W2 = 1 as the reference transmitted 
message. Let Ax^q(w,) = x^q(w/) -x"q(1), as lO, with h replaced by (/zi,/z2), for 1 <m<M, 
we have the following results 

P{m = m) <p{{{huh2) i 0^'-i},{y™o"i i Un-i} 

+p{{{hM)eO'^-'},{{h,M)i(yi} 

where the definition = (|), = (|) are again used. 
For 1 < m < M, given a {h\,h2) and a codebook, 

(49) 

+ £ /4|/7iAx'[;o(wi)+/i2Ax^oMl'<4mr(l +5)0,2} 

vv/(l,l) 

where Id{) is the indicator function. Different from [fT4j. we do not average over the random 
coding ensemble here from (1491 ). Instead, we use the following inequality 

Id {\hx^•s!^^Q{wl)+h2^xlQ{w2)\^ < 4mr(i + 6)o2) 

1 [ |/ziAx'/^o(wi)+/z2Ax^"o(w2)n (50) 
^"^^^'Py \mT{l + h)ol 

Then we turn to upper bound the term P{E\Er^va = m) 



(48) 



B. Analysis of P{E\Er,m = m) 

We consider the GLRT decoder at the destination. This decoder has no knowledge of decision 
time m. Thus the receiver has to decode both the decision time and source information message. 
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The P{E\Er,i!n = m) is upper bounded by 

P{E\Er,m = m) < P{{{gi,g2,gr) G O^}) 

+ P{E, {{gug2.gr) i Ol}\Er,m = m) 

and similar to [fT4l . we have 

/'({(l,l)^M)}|£„m = m) 

M _ 

< E P{p{y%\l.m)<p{y'f^^\w,m')\Er,m = m) 



(51) 



(52) 



m' = 1 



and for m' > m, (the case for m' < m can be derived by interchanging the role of m' and m, thus 
omitted here ), the term inside the sum (|52] ) can be upper bounded by 



(53) 



^iAx-o(wi)+^2Ax-o(#2) 

glAxfAm)+g2AxlJw2) +grAx^^,iw) 



C. Averaged over the Gaussian random ensemble 

We are ready to average the term P{m = m)P{E\Er,'m = m) over the Gaussian random 
ensemble. From (|48]).(|5n). for 1 < m < M, it can be upper bounded by sum of the four terms 
discussed below, we will show that the four terms have exponents larger than or equal to 

dm,D{r) +d„r^R{r) : 

(1) Pi{{gug2.gr) e 0"^})p(^{ihuh2) e 0'^-'},{ihM) i O'^}): 

It is clear this term has the exponent of d,„,i){r) +dm,R{r)- 

(2) P{E,{{gi,g2,gr) i 0"^}\Er,m = m)p(^{{hM) i O^-^}, 

{y'^Q^^ ^ lim-ijj: This term is more involved. For the case m = 1, the term {2) becomes to 
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zero. For 1 < m<M, given hi, gi, from (|48] ). (|50] ). (152] ). (|53T ). we have 
p({yTo'^nn-i})-P{{{lA)^w}\-Er,m = m) 



< 



.^7^(1,1) 



+(1+5) 



(ra-l)r -(ra-l)r5 



|/ziAx';^oHwi)+/i2Ax^;onv^2)| 
(4(m-l)r(l+6)a2) 



|/^iAx7oi(wi)+/r2Ax^;Q^(w2)| 
(4(m-l)r(l+5)a2) 



lvlV(4a2 



m— 1 ; 



w^(l,l) 
M 

e 

m'=l 



£ .-|vlV(4o2) 



|/iiAx7,o1(wi)+/i2Ax^oHm^2)| 



m— 1 / 



w/(l,l) 
M 



m'=l 



(4(m-l)r(l+5)o2) 
. £ g-|z„,,|V(4(m-l)r(l+6)a2) 

|/iiAx';^oi(wi)+/z2Ax^^oi(w2)f 
(4(m-l)r(l+5)a2) 

M 

£ g-|z,„,|V(4(™-l)r(l+8)<^?) + £ g-|z;,^,|V(4(m-l)r(l+6)o?; 



= 1^ ^'■^p 

W7^{(l,l),w} 
M 



m'=l 



where 



JT A 



(54) 



(55) 



(56) 



(57) 



(58) 



§;^(/?lAx'[;ol(M)i)+/r2Ax^,o'(^2)) Zm' 

Note by separating the term of m> from the first sum in (|58T ). then the product of the first two 
sum in (|58l) can now be averaged separately. We first consider the single term 

„,|V(4('n-l)r(l+5)a?) 



(59) 



Let 



and 



X/ = [A^l,,(wi),A^2,Kl^2),^r,Kl, l),A^r,Kl^)]^, I < i < MT 



X — [Xf ,X2 , ...X^T']^ 



It can be verified that |z^p = X^RX, where R is a block diagonal matrix of the form R 
diag{Ri,...RMT)- 
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For 1 < z < (m- l)r, 



Ri = C"C 



where h'; = %hi and 





h'2 








gi 


gi 

































For (m-l)r+l <i<mT 



For niT + 1 < z < m'T 



For m'T + \<i< MT 



R/ = ga*g«, ga = [^l,,?2,0,0] 



R/ = gf*gfc, gZ. = [.?l,,?2,.?r,0] 



R/ = gf * gc , gc = [gl ,g2,0, gr 



Let W = diag(R\, ...,R^j^), D = diag(Di, ...,Dmt) where denotes R,- with the 5-th row and 
the 5-th column replaced by zero vector. (when s = 0, nothing is changed). D,- = diag{2P, IP, P, IP) . 
is similarly defined. Averaging (|59] ) over the Gaussian random ensemble, we have 

1 



det 1 + 



(4(m-l)r(l+5)o2) 



(60) 



Let k 



((m-i)r(i+6))' 



we have 
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(/ j For {wi ^ 1,#2^ 1}, * = , 

1 



det(l+«^) dct(l I ^pc/cy^'^^^ 



1 



1 + 



fcp(|glP + |g2p 



1 



1 + 



fcp(2|gip+2|g2p+|g,p) 



1 



(m'—m)T 



1 + 



fcp(|glP + |g2p+|g,.p) 



{M-m')T 



<- 



[l+Pi\hi\^ + \h2\^+\gl\^ + \g2 



|2N](m-l)r 



[l+P(|.?lP + l.?2 



[l+P(klP + k2|2+|g.P)](^-'")^ 

Notice that (|6T1) does not depend on m' and follows from the inequality 

l + ^i\hi\^ + \h2\^+\gi\^ + \g2\' 



kpC"C 



det I + 



> 



(61) 



Define \hi\^ = p""', \gi\^ = p^P'. Use the union bound for 'B„^, summing over all m' = 1, ...,M and 
over all messages {vvi ^ l,W2 ^ I}, and average over the channel realizations where {{hi,h2) ^ 
{{gi,g2,gr) ^ 0'£}, use the techniques developed in |l3l, we obtain the exponent of this 
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correlated term as 



d'=Lir) = inf {ai + a2 + Pi + p2 + p. + ry;r°(a„ p„ r) } 

> inf {ai + a2 + pi+p2 + P,+r/r°(a/,P/,r)} (62) 

a/,p;>0, 

> inf {ai + a2 + r/r°(a„p„r)} (63) 

C XnPi>0 , 

+ inf {pi+p2 + P. + r/r'(a,-,p„r)} (64) 
c x,-,[5,>0 , 

>dm,R{r)+dmn{r) (65) 



where 



(66) 



= (m-l)max[(l-ai)+,(l-a2)+,(l-pi)+,(l-p2)+] 

+ max[(l-pi)+,(l-p2)+] 

+ (M- m)max [(1 - pi)+, (1 - ^3) + , (1 - Pr)+] - rM 

(l63l) . (|64|) follows form the fact that the infimum of diD occurs where /,^=*^(a/, P,, r) if T 
is large enough(i.e. T > 4). Notice the event 



{(a,) : (m-l)max[(l-ai)+,(l-a2)+]-rM>0} e O"^-^ (67) 



It can be checked that ([63]) > dm.R{r), and ([64]) > dm,D{r)- 
(II) For {wi^l,W2 = l},s = 2 , 



1 



det(H-^) [l + p(|/.iP+|^i|2)](-i)- 

1 



[1 + Pk 



1 



(68) 
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Apply similar arguments, we obtain the exponent of this correlated term as 
C.wl'-) = inf {ai + pi + p. + r/r'(a„ p„ r) ] 

> inf {ai+pi+p, + ry^r'(0CoPn'-)} (69) 

C X/,P,>0 , 

an-' 

> inf {ai+r/,r'(a„P„r)} (70) 

a/,p,>0, 

+ inf {pi+p, + r/^=2(a,,p,r)} (71) 

«/,P,>0 , 

> dm,R{r) + d,n,D{r) (72) 

/,r'(a„p„r) = (m-l)max[(l-ai)+,(l-pi)+] 

+ max(l-pi)+ (73) 
+ (M - m) max [(1 - pi)+, (1 - pr)+] - 

|ai : (m-l)max(l-ai)+-^M>o| G O^'"^ (74) 

The case for {wi = l,'W2 ^ I}, s = I is similar and omitted here. Thus this correlated term can 
be ignored with respect to the term (i) in DMT analysis. 
For averaging the first sum in (|58T ). define 



where 



Note 



h'-^h' 1 <?■< (m-l)r 
otherwise 



where 



h[ h'2 



h' 

and s is chosen according to w as before, we have 



Averaged by {{huhj) ^ O'^^^}, {(^i,^2,^r) ^ O^}, the exponent c?„,(/,,.)(r) can be derived. 
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To average the second sum in (|58] ). let Q'^ = W, except for 1 < z < (m — l)r, Q- = R^j, use 
the union bound over w, by averaging the term 



1 



w^(r,i)det(l + ^ 



(76) 



over {(/?i,/z2) ^ }, {{gi,82:gr) ^ O'^}^ we obtain the exponent J,„(g.)(r). Analogous to 
Appendix IA-A[ Appendix lA-Bl it can be shown that d^ (j^.^{r)>dm.R{r), d^ (^g.){r) = dm.pi^)- 
Collecting all the results above, the term (2) can be ignored compared to the term (ij in DMT 
analysis. 

(3) P{{{gug2,gr) e 0'^})p{{{hiM) i 0^'"U,{y"o"' i U^-y}): These two terms can be av- 
eraged over the Gaussian random ensemble separately. It is clear the first term has the exponent 
of dm.pir), and the second has d^j^ij.^{r). 

(4) P{E,{{gug2,gr) i OlWr.v(y = m)p{{{hyM) e 0^-1} , 

{(^1,^2) ^ O^})'- These two terms can be averaged over the Gaussian random ensemble 
separately. It is clear the first term has the exponent of (g.)(r), and the second has dm,R{r). 

Appendix C 
Proof of Theorem [3] 

To upper bound the second term on the right hand side of (1261 ). we further partition the error 
event into two error events depending on the MAF's outage event Omaf of a,, (3,, where the 
MAP is being used, 



M M 

P [ E.Omaf, < {{hiM) i 01 }, {y,?^o iUu})]^ (77) 



P \E,OMAF.\{{hxM) i 01 },{y,^o i Um} 

Notice we do not apply the standard DMT analysis [9] here by simply using P{Omaf) to 
upper bound the first term on the right hand side of dTTt (denoted as P^^\^ for notation ease), 
otherwise the lower bound of DMT we obtain for (|77l ) will be the same as MAF protocol, only 
2-^, not 2-r. 



25 



For the the analysis of P^^' , note the fact 



we have 



P{ E 



p'''ma,,%)eOMAF) 



a/,|3,-) e Omaf. { {{hM i Oj},{yJo i Um} 



< 1 



<P 



(a/,|3,) G Omaf 



1 



{{huh2)iO^],{y^^^iUM] 
Recall (|49l) . (I50l) . by averaging over the Gaussian random codebook, we have 

f'"(5|(ai,|5,)eO„^f) 



<£p-*?i(.-.,)-;i + p-^[(.-»i.(«,..,)-: 

i=\ 



(78) 



Finally, averaging the right hand side of (|78T ) over (a/, (3/) G < Omaf fl ( (0CI7OC2) ^ C)^^ ) k foi" 



T is large enough( but finite), the derivation of = p ^ becomes identical to (|29l) . (|3T| ) 

and dpist{r) > 2 — r. Then we tum to upper bound the second term on the right hand side of (fTTl) . 
denoted as p2«rf(^). As (ES]), p2nrf(^) < i:/P2«rf(^/). Use Pmaf (£,Omaf), where MAF is always 
used, to upper bound /'^"'^(^), we have ^'^"'^(^(1,2)) > 3 — r, however, only /'^"'^(^/) > 2 — y is 
obtained, / = 1,2. To get a tighter upper bound for /'^'"'(^i), (|49l) , (|50l) are again used and apply 
the techniques developed in Appendix IB-CI to MAF protocol, where the equivalent MAF model 
can be expressed as 



xi+ n + v (79) 

'2h\b\MT/2 g\^MT/2\ [^2^^7/2 

and b is chosen to be of exponential order zero [|6l, [l3l, which satisfies the relay's transmission 
power constraint. We have 

Pi2"^({(l,l)^M)})<e-l^'l'+ 



M 



M 



52 exp 

w^{(l,l),vv} 



|/ziAx2o(wi) +/z2Ax2=o(w2)h 



(80) 



(2Mr(l+5)a2) 

where P\"^ ({(1, 1) — t- vv}) denotes the probability that (1,1) is decoded as w for type-1 error in 
the event of the second term in (|77]) given /z,,^, and z and z' are redefined in this section as 

£lMr/2 



Axi(vv) 



(81) 
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i! 



^1 T , n 



Axi(w) (82) 



Apply the union bound and average them over Gaussian random ensemble and corresponding 
channel realizations. For the product term on the right hand side of (|80l) . it can be averaged 
separately over the Gaussian random ensemble, resulting in the exponent dproduct{f) identical to 
(|29l ). For the first term on the right hand side of (|80l l. we obtain its exponent equal to 

f MT 1 

inf <^ai+a2 + |3i + (32 + |3,+ — fc(a,-,|3,-,r) ^ >2-r (83) 

where 

fc(a,-,p,-,r)=max{2(l-pi),l-(|3, + ai),2-(ai+pi)}-r 



and C= joMAFfl (^(oci, 0(2) ^ | 
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